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ABSTRACT 


Modeling  and  simulation  have  proven  invaluable  to  the  design  of  various  systems 
and  planning  of  many  experiments  since  the  free-electron  laser  was  invented  in  the 
1970s.  This  thesis  illustrates  methods  used  for  these  tasks  and  for  the  validation  and 
development  of  further  theory.  PEL  systems  and  their  components  are  described,  briefly 
discussing  how  they  fit  into  design  considerations.  Basic  theory  of  PEL  operation  is 
reviewed,  including  the  resonance  condition,  PEL  wave  and  pendulum  equations,  and  the 
concepts  of  undulator  taper,  extraction  and  beam  quality.  Two  computer  programs  used 
for  PEL  simulation  {FEL3D/FEL4D,  developed  at  the  Naval  Postgraduate  School,  and 
Ginger,  developed  at  Lawrence  Berkeley  National  Laboratory)  are  discussed  and 
compared,  then  applied  to  specific  PEL  cases.  The  cases  studied  serve  to  evaluate  key 
performance  characteristics  of  PELs.  The  final  case  study  directly  compares  output  from 
the  two  programs. 

In  addition  to  this  type  of  physics  simulation,  computer  tools  are  also  vital  for 
component  design.  Relevant  methods  are  illustrated  using  the  example  of  an  electrostatic 
cathode  test  cell  and  drawings  made  using  the  Los  Alamos  National  Laboratories  design 
tool,  Poisson  Superfish.  After  reviewing  the  principles  involved  and  system  constraints, 
a  new  test  cell  is  designed  for  use  with  future  NFS  experiments.  The  effects  of  anode- 
cathode  spacing  and  anode  aperture  size  are  examined  for  two  different  anode 
configurations  and  recommendations  are  made  for  future  design  iterations. 
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I. 


INTRODUCTION 


Free-electron  lasers  (FELs)  are,  like  any  other  lasers,  sources  of  coherent 
electromagnetic  radiation  (light),  which  is  amplified  by  stimulated  emission  from 
electrons.  Free-electron  lasers  differ  from  other  types,  however,  in  how  the  light  is 
generated.  Rather  than  relying  on  electrons  bound  within  matter,  FELs  use  a  high-energy 
beam  of  electrons  in  vacuum.  They  are  distinguished  from  solid-state  and  gas  lasers  by 
the  flexibility  of  their  designs:  Adjustment  of  a  few  system  parameters  completely  alters 
their  output  wavelength,  whereas  conventional  lasers  are  restricted  by  the  quantum 
properties  of  their  lasing  media.  Some  of  the  relevant  FEE  parameters  will  be  set  during 
system  construction;  but  others,  such  as  the  electron  beam  energy,  can  be  altered 
continuously,  allowing  wide  flexibility  in  application. 

Years  of  study  with  FELs  have  also  shown  the  systems  to  be  powerful,  reliable 
(98%  or  more  uptime  on  many  systems)  and  efficient  (10%-20%  wall-plug  efficiency 
according  to  some  estimates).  On  the  other  hand,  FEE  systems  have  historically  been 
large  and  expensive;  a  smaller  footprint  and  the  ability  to  be  supported  on  mobile 
platforms  are  active  areas  of  research  for  some  applications. 

Modeling  and  simulation  have  proven  invaluable  to  the  design  of  various  systems 
and  planning  of  many  experiments  since  the  free-electron  laser  was  invented  in  the 
1970s.  This  thesis  illustrates  methods  used  for  these  tasks  and  for  the  validation  and 
development  of  further  theory. 

Chapter  II  describes  FEL  systems  and  their  components,  briefly  discussing  how 
they  fit  into  design  considerations.  Chapter  III  reviews  the  basic  theory  of  FEL 
operation,  including  the  resonance  condition,  FEL  wave  and  pendulum  equations,  and  the 
concepts  of  undulator  taper,  extraction  and  beam  quality.  Chapter  IV  discusses  and 
compares  two  computer  programs  used  for  FEL  simulation  {FEL3DIFEL4D,  developed 
at  the  Naval  Postgraduate  School,  and  Ginger,  developed  at  Lawrence  Berkeley  National 
Laboratory)  and  applies  them  to  study  specific  FEL  cases.  The  cases  studied  serve  to 
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familiarize  readers  with  the  simulations’  output  and  highlight  key  performance 
characteristics  of  FELs.  The  final  case  study  directly  compares  output  from  the  two 
programs. 

Chapter  V  delves  into  the  specifics  of  physical  design,  using  the  example  of  an 
electrostatic  cathode  test  cell  and  drawings  made  using  the  Los  Alamos  National 
Laboratories  design  tool,  Poisson  Superfish.  The  principles  and  design  constraints  for  a 
test  cell  to  be  used  with  future  NFS  experiments  are  examined.  The  effects  of  anode- 
cathode  spacing  and  anode  aperture  size  are  studied  for  two  different  anode 
configurations  and  recommendations  are  made  for  future  design  iterations. 
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II.  FREE-ELECTRON  LASER  SYSTEM  DESCRIPTION 


Free-electron  lasers  can  have  two  basic  configurations:  oscillator  or  amplifier. 
Most  components  are  common  to  both  designs,  as  shown  in  Figure  1.  The  difference  is 
how  the  light  is  extracted  from  the  electron  beam — amplifiers  amplify  an  existing  optical 
field  over  one  pass  while  oscillators  use  optical  resonators  similar  to  conventional  lasers. 
The  following  sections  discuss  the  purpose  and  general  properties  of  the  system 
components.  Figure  1  depicts  a  particular  FEL  system  with  a  few  illustrative  properties: 

1.  Separate  beam  paths  for  running  in  oscillator  or  amplifier  mode.  This 
shows  the  difference  between  the  two,  but  the  dual  beamline  configuration  can  be  found 
in  real  systems,  such  as  the  experimental  setup  at  Thomas  Jefferson  National  Laboratory. 

2.  A  superconducting  radio-frequency  linear  accelerator,  and  the  helium 
refrigeration  plant  associated  with  it.  Normal-conducting  systems  are  also  used  widely. 
Electrostatic  (DC)  and  storage  ring  accelerators  have  also  been  used  for  FELs. 

3.  Beam  recirculation  for  energy  recovery  (see  “Beam  Disposal”  below). 


Superconducting  Etectron  Accelerator 
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Figure  1 .  FEL  System  Schematic 
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A.  INJECTOR 

The  injector  generates  the  beam  of  free  electrons  from  a  cathode.  Cathodes  come 
in  several  varieties  and  generally  produce  electrons  either  via  thermionic  emission,  the 
photoelectric  effect  or  field  emission.  Thermionic  cathodes  emit  electrons  when 
electrically  heated  to  a  point  where  the  thermal  energy  in  the  metal  exceeds  its  work 
function.  Photocathodes  emit  electrons  due  to  the  photoelectric  effect.  When  they  are 
illuminated  by  a  drive  laser,  as  long  as  the  energy  per  photon  exceeds  the  work  function, 
electrons  are  liberated.  Field-emission  cathodes  are  not  necessarily  metal,  and  use 
geometry  to  create  points  where  high  electric  fields  exist  once  an  external  potential  is 
applied.  For  example,  diamond  field-emitter  array  (DFEA)  cathodes  are  wafers  of 
diamond  film  grown  on  a  suitable  substrate  such  that  the  diamond  forms  an  array  of 
microscopic  pyramids;  the  field  is  enhanced  enough  at  the  pyramids’  tips  that  electrons 
escape  by  Folwer-Nordheim  tunneling. 

Regardless  of  cathode  type,  injectors  will  apply  some  external  field  to  the  emitted 
electrons,  accelerating  them  to  modest  energies  before  they  enter  the  accelerator 
component.  There  will  also  be  some  amount  of  instrumentation  and  beam  transport 
apparatus  (e.g.,  focusing  magnets)  to  prepare  the  beam  for  entry  to  the  accelerator  and  to 
perform  initial  characterization  of  the  electron  beam. 

Injector  design  can  be  important  to  overall  system  performance  in  that  it  must 
provide  electrons  to  the  accelerator  with  optimal  bunch  length  and  repetition  frequency. 
Furthermore,  the  beam  quality  (i.e.,  the  spatial  or  angular  spread  of  the  electrons  or  the 
relative  disparity  of  their  kinetic  energies;  see  discussion  in  Chapter  III)  provided  by  the 
injector  is  the  highest  that  will  be  present  in  the  entire  system. 

B.  ACCELERATOR 

The  electron  beam  is  next  passed  through  an  accelerator,  which  raises  its  energy 
to  the  desired  level.  Direct-current  (DC)  acceleration  is  possible,  but  high-energy 
applications  usually  use  radio-frequency  (RF)  fields.  Electron  accelerators  generally 
consist  of  a  series  of  evacuated  metal  cavities  to  which  external  RF  power  is  applied  in 
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order  to  generate  strong  eleetromagnetie  fields  inside.  The  alternating  fields  ean  be  timed 
so  as  to  provide  net  positive  acceleration  to  the  sequence  of  electron  pulses  provided  by 
the  injector. 

C.  UNDULATOR 

The  undulator  consists  of  a  series  of  alternating  transverse  magnetic  fields,  which 
produce  transverse  accelerations  of  the  electron  beam.  This  causes  the  electrons  to 
radiate,  producing  spontaneous  emission  like  population  inversion  does  in  a  conventional 
laser.  The  transverse  motion  of  the  electrons  also  allows  for  transfer  of  energy  between 
the  electrons  and  the  optical  field,  stimulating  further  emission,  creating  a  coherent  light 
source.  The  relativistic  speeds  of  the  electrons  ensure  this  optical  field  takes  the  form  of 
a  narrow  coaxial  beam. 

Undulators  can  be  linear  or  helical  in  their  polarization.  Linear  undulators  (see 
Figure  2)  consist  of  a  series  of  dipole  magnets,  alternating  in  orientation,  which  create  a 
magnetic  a  field  that  is  periodic  in  one  transverse  plane. 


Figure  2.  Linear  Undulator,  showing  (1)  Magnets,  (2)  Electron  Beam  and  (3) 

Synchrotron  Radiation  (not  discussed)  [After  1] 

Helical  undulators,  as  shown  in  Figure  3,  usually  consist  of  permanent  magnets  or 
electromagnets  coiled  around  the  beam  path  in  a  helix,  and  provide  a  magnetic  field  that 
is  periodic  in  both  transverse  directions  (see  Chapter  III).  There  are  also  planar-geometry 
undulators  that  can  provide  a  helical  magnetic  field,  such  as  the  APPLE  (Advanced 
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Planar  Polarized  Light  Emitter)  system  [2],  In  the  APPLE,  two  linear  undulators  are 
operated  with  their  planes  of  polarization  at  an  angle  to  one  another,  thus  producing  net 
polarization  that  can  be  helical  (elliptical  or  circular)  or  linear  in  a  specified  direction, 
depending  to  the  angle  used. 


Z 


Lo)N-Oarbon  Sfee/  Pole 


OD-B.3mm 


Figure  3.  Helical  Undulator  [From  3] 


D.  OPTICAL  COMPONENTS 

In  FEL  oscillator  designs,  the  undulator  is  contained  in  a  resonator  cavity,  with 
one  fully  reflective  and  one  partially  transmissive  mirror.  As  the  light  repeatedly  transits 
this  resonator,  it  serves  to  stimulate  further  emission  from  new  electrons.  On  each  pass,  a 
portion  of  the  light  is  outcoupled  via  the  partially  transmissive  mirror.  In  most  FEL 
amplifier  designs,  the  initial  optical  field  is  produced  by  a  seed  laser  that  is  amplified 
while  passing  only  once  through  the  undulator.  There  are  also  self-amplifying 
stimulated-emission  (SASE)  designs  that  use  the  spontaneously-emitted  light  from  the 
first  portion  of  the  undulator  as  a  seed.  Amplifier  designs,  in  general,  require  longer 
undulators  than  oscillators,  with  SASEs  requiring  significantly  more  length. 

E.  BEAM  TRANSPORT  AND  DIAGNOSTICS 

The  physical  layout  of  the  FEL  system  will  determine  exactly  which  components, 
such  as  bending  magnets  and  focusing  magnets,  are  required  to  properly  transport  the 


6 


beam.  Most  systems  will  also  have  several  types  of  diagnostics  at  various  places  along 
the  beamline,  both  for  application-specific  purposes  and  to  track  the  results  of  beam 
transport. 

F.  BEAM  DISPOSAL 

After  light  is  generated  in  the  undulator,  one  must  dispose  of  the  spent  electrons. 
Ultimately,  the  beam’s  energy  can  be  dissipated  into  something  as  simple  as  a  metal 
target  as  long  as  appropriate  cooling  and  radiation  shielding  are  provided.  This  is  labeled 
as  “Beam  Dump”  in  Figure  1.  Since  light  generation  only  depletes  a  small  fraction  of  the 
electrons  beam’s  energy,  some  designs  recirculate  the  beam  through  the  accelerator;  if  re¬ 
injected  with  the  proper  phase,  the  electrons  will  give  up  significant  energy  to  the  RF 
fields,  boosting  system  efficiency  and  making  beam  disposal  much  less  hazardous.  This 
technique  is  referred  to  as  “energy  recovery,”  and  linear  accelerators  that  use  it  (not  all  of 
which  are  FELs)  are  known  as  “energy-recovery  linacs”  or  ERLs. 
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III.  BASIC  FREE  ELECTRON  LASER  THEORY 


The  centerpiece  of  PEL  physics  is  the  interaction  between  electrons  and  light  that 
takes  place  within  the  undulator.  The  physics  of  particle  acceleration  and  laser  beam 
transport  are  vital  to  the  operation  of  the  greater  system;  but  the  undulator  is  where  the 
interaction  characteristics  of  ffee-electron  lasers  are  determined. 

It  is  helpful  to  begin  a  mathematical  description  of  these  interactions  by  defining  a 
few  useful  coordinates: 

_  ct 
^  ~  L 

=  (k  +  k^z  — {3.1} 

y^^  =  L\(k^k,)P^-k\ 
dz 

The  normalized  time  z  serves  as  a  temporal  variable  along  the  longitudinal  spatial 
coordinate;  z  =  0  corresponds  to  the  beginning  of  the  undulator  and  r  =  1  to  its  end.  The 
electron  phase  ^  measures  an  electron’s  relative  position  in  a  portion  of  the  electron 
beam  the  length  of  one  optical  wavelength  A,  and  the  electron  phase  velocity  v 
measures  its  rate  of  change.  Here  z  is  the  electron’s  position  along  the  undulator; 
CO  =  kc  is  the  frequency  of  the  optical  field;  L  is  the  length  of  the  undulator; 
k  =  2n !  A  and  k^=lK  t  A^  are  the  wave  numbers  corresponding  to  the  optical 

wavelength  A  and  undulator  period  A^  respectively;  and  P^=v^ !  c  is  the  normalized 

longitudinal  electron  velocity.  We  will  see  that  these  definitions  play  an  important  role  in 
the  free-electron  laser’s  interaction. 

A.  RESONANCE  CONDITION  AND  PENDULUM  EQUATION 

The  electrons’  motion  is  determined  by  the  Lorentz  force  equations,  namely 

X  5) ,  and 

dt  mc^  '  2} 

dt  me 


9 


For  a  helical  undulator  (see  Chapter  II),  the  magnetic  field  can  be  expressed  as 

Knd  =  B{cos[k^z],sm[k^z],0),  {3.3} 

while  the  optical  fields  will  take  the  form 

E  =  £'(cos  y/,  -  sin  y/,  0) 

=  £'(sin  y/,  cos  y/,  0),  where  {3.4} 

y/  =  kz  -  cot  (p. 

In  the  preceding  equations,  E  is  the  optical  electric  field  strength,  B  is  the  undulator’ s 
magnetic  field  strength,  cp  is  the  field  phase,  and  e  and  m  are  the  charge  magnitude  and 
mass  of  an  electron.  In  Equation  {3.2},  the  pertinent  magnetic  field  is  the  sum 

^  =  ^und  +  ^opt  ■ 


Substituting  {3.3}  and  {3.4}  into  {3.2}  and  assuming  that  /?^«1  gives  the 
transverse  motion  as 


= - (cos[A:oz],sin[A:oz],0),  where 

r 


iTimc 


2  ■ 


{3.5} 


The  “undulator  parameter”  .ST  is  a  dimensionless  measure  of  field  strength  and  is  of  order 
unity.  Note  that  this  derivation  assumes  ideal  injection  into  a  helical  undulator  such  that, 
when  solving  for  position  as  a  function  of  time,  the  constants  of  integration  are  all  zero. 
In  the  helical  case,  B^^  =  B ;  but  a  similar  analysis  holds  for  linear  undulator  geometry 

with  B^^=B/42. 


The  evolution  of  electron  energy  can  be  analyzed  through  changes  in  y  and  in  v. 
In  particular,  from  the  second  equation  in  {3.2}  and  Equation  {3.5}: 


dr  eKE  , 

—  = - cos(^  +  (p). 

at  ymc 


{3.6} 


The  electron’s  velocity  is  related  to  its  energy  and  hence  its  Lorentz  factor  y,  specifically, 


Y-^  =  \-Pi-Pt. 


{3.7} 
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Equation  {3.5}  shows  that  the  magnitude  of  the  transverse  velocity  is 
merely  /  y  ;  thus,  for  y  »  1 ,  Equation  {3.7}  reduces  to 


«1- 


\  +  K^ 


{3.8} 


As  an  aside,  we  consider  next  the  condition  under  which  optimum  energy  transfer 
occurs  between  the  electron  beam  and  optical  field — ^referred  to  as  “resonance.”  In  this 
state,  each  electron  emits  (or  amplifies)  one  wavelength’s  worth  of  light  (A)  for  each 
undulator  period  (distance  A^)  through  which  it  passes.  If  the  interaction  between  the 
electron  and  the  optical  field  is  thought  of  as  a  race  down  the  undulator,  then  the  light 
leads  in  velocity  by  c(l  -  /?^)  and  wins  by  a  distance  of  A  in  one  undulator  period.  The 
time  passed  is 

A/  =  A.  {3.9} 

cP, 


The  winning  distance  can  be  expressed  as  the  velocity  difference  multiplied  by  the  race 
time: 

1-A 


A  =  c{\-  =  /Ifl- 


A 


{3.10} 


From  Equation  {3.8},  we  see  that  the  condition  for  resonance  is 


A 


P-. 


,  or,  for  y  »  1, 


A  ~  A, 


i  +  t:" 


{3.11} 


which  relates  the  FEL’s  output  wavelength  to  fundamental  properties  of  the  electron 
beam  and  undulator.  This  important  relation  highlights  the  tunable  and  designable  nature 
of  FELs.  Adjustment  of  any  or  all  of  the  system  parameters  that  impact  K  allows 
designers  to  build  systems  that  emit  at  any  of  a  wide  range  of  wavelengths,  or  even  to 
adjust  K  on  an  existing  system. 

Using  the  “resonance  condition,”  Equation  {3.11},  in  the  definition  of  phase 
velocity,  and  recognizing  that  k  ::$>  (i.e.,  y  »  1 )  shows  that 
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V  =  =  L 

dr 

ItiL 


dv 


A 


(k  +  k^) 
l  +  K 


V 


-k 


j 


(-2)^  = 

Y 


\  K  j 


'^-Lk 

dy 

r  ' 


y 


,  and 


{3.12} 


From  the  definition  of  undulator  period  (L  =  N  where  N  is  the  number  of  periods), 
we  find 

dv  =  A7rN^.  {3.13} 

r 


Differentiating  v  with  respect  to  r  and  using  the  result  of  Equation  {3.6}  gives: 


dv 

dr 

\a\ : 


d^C 

=  |fl|cos((^  +  (p),  where 


AKNeKLE 

2  2  ■ 

Y  me 


{3.14} 


This  is  known  as  the  “FEE  pendulum  equation”  for  its  similarities  to  the  equation  of 
motion  for  a  simple  pendulum.  Using  this  equation,  the  evolution  of  FEE  electrons  is 
often  described  with  phase  space  {<^  vs.  v)  graphs. 

Parameter  a  is  a  dimensionless  measure  of  the  laser  field  amplitude  E.  Values  of 
a  <  TT  are  considered  the  “weak  field”  regime,  while  a  »  .tt  is  considered  “strong 
fields.”  The  FEE  theory  as  presented  is  valid  for  both  regimes,  as  well  as  intermediate 
values. 


B.  FEE  WAVE  EQUATION 


The  fee’s  optical  field  will  evolve  according  to  Maxwell’s  wave  equation: 


V" 


V 


JL  ^ 

A{x,t)  = - 

c 


{3.15} 


In  this  equation,  A  is  the  optical  vector  potential  and  J  is  the  transverse  current  density 
due  to  the  undulations  of  the  electrons.  The  vector  potential  is  related  to  the  optical  fields 
by: 


E  = - , 

c  dt 


B  =  VxA. 


{3.16} 
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{3.17} 


Assuming  a  complex  optical  electric  field  in  the  form 

where  a  =  kz- cot  represents  the  earrier  wave,  the  veetor  potential  becomes 

A{x,t)  =  ^^^nE\  {3.18} 

k 

where  n  is  the  polarization  vector.  Given  that  the  light  speetrum  produeed  is  not  only 
finite  but  narrow  (i.e.,  a  laser),  it  is  safe  to  assume  that  the  amplitude  and  phase  vary  only 
slowly  along  the  seale  of  an  optical  wavelength  {A  =  27r !  k)  and  the  scale  of  the  optieal 
frequeney  {co  =  kc).  In  this  case,  one  can  assert 

^^(o\e\,  ^<^(0(p,  ®«A:|.E|,  ^^kcp.  {3.19} 

dt  dt  dz  dz 

Substituting  Equation  {3.18}  into  Equation  {3.15}  and  using  conditions  {3.19}  to 
neglect  terms  with  two  derivatives,  the  relation  simplifies  to 

Vl+2ik(—  +  -—]  E{x,t)  =  -^^(j  {3.20} 

\dz  c  dt  Jj  c  ^  ' 

where  the  transverse  Eaplaeian  symbol  (V^)  denotes  differentiation  only  in  the  x  andy 

coordinates.  In  this  new  parabolie  form  of  the  wave  equation,  the  souree  term  can  be 
expressed  as 

J{x,t)  =  [x-r.{t)].  {3.21} 

i 

In  this  expression,  r.  (t)  is  the  position  of  the  /th  eleetron  at  time  t  and  is  the 

transverse  veloeity  diseussed  in  Equation  {3.5}.  For  a  helical  undulator  producing 
circularly  polarized  light,  n  =  (-/,  1,0) ,  thus 

—  ^  K  ^ 

=  Re - .  {3.22} 

\  r  J 

Invoking  the  method  of  eharacteristics  to  introduce  the  coordinate  u  =  z-ct 
(whieh  follows  the  light  through  the  evolution),  the  derivatives  in  Equation  {3.20}  are 
transformed  as  follows: 
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{3.23} 


d  _  d 
dz  du 
d  _  d  d 

dt  du  dt 


Recalling  r  as  defined  in  Equation  (3.1),  Equation  {3.23}  becomes 


V=  +  .2‘*  ® 


L  dr 


E{x,  t)  =  -AnieKk 


-ii 


\  r 


)p{x,t). 


{3.24} 


where  the  notation  <expression>  indicates  an  average  of  the  expression  over  all  electrons 
and  the  electron  particle  beam  density  in  a  small  volume  element  dV  is 

p{x,t)  =  \Y,S^^^[x-fXt)]dV.  {3.25} 

V  i 


To  simplify  presentation  (and  mathematics),  we  utilize  the  definition  of 
dimensionless  laser  field  amplitude  from  Equation  {3.14}  and  introduce  the 
dimensionless  current  density 


%Np{KeKLf 


{3.26} 


Values  of  j  <k  are  “low  gain”  cases,  while  7  »  is  the  “high  gain”  regime;  as  with 
the  values  of  a,  this  FEE  theory  is  valid  for  all  values  of  j. 

Transforming  transverse  coordinates  to  dimensionless  form  by  x  =  x^Jk  /  (2L) 
and  y  =  y^jk  /  (2L) ,  the  final  form  of  the  FEE  wave  equation  is  obtained: 

i(x,0  =  -(je“'^).  {3.27} 


-V\+  — 
4  dr 


The  modified  Laplacian  symbol  is  used  here  to  denote  differentiation  with  respect  to 
the  normalized  coordinates  x  and  y .  The  terms  of  Equation  {3.27}  involving  this 
spatial  derivative  describe  the  evolution  of  the  optical  field  due  to  diffraction  effects. 
Terms  involving  the  time  derivative  describe  its  evolution  due  to  gain  from  the  electron 
beam.  The  term  <  je~'^  >  is  a  measure  of  how  strongly  bunched  the  electrons  are, 
showing  that  this  portion  of  the  dynamics  have  a  significant  effect  on  optical  output.  It  is 
worth  noting  also  that  in  Equation  {3.27},  the  evolution  of  a  depends  on  electron 
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distribution  ^ ,  while  in  {3.14}  the  evolution  of  ^  depends  on  |fl|;  thus,  the  FEL  wave 
and  pendulum  equations  form  a  feedback  loop.  Bunched  electrons  amplify  light  better, 
while  the  optical  field’s  coherent  evolution  promotes  electron  bunching.  The  gain  from 
this  process  only  grows  for  so  long,  however,  as  electrons  will  eventually  overbunch,  at 
which  point  the  FEL  reaches  saturation. 


C.  UNDULATOR  TAPER 


In  high-gain  FELs,  saturation  can  be  reached  before  the  end  of  the  undulator, 
lowering  overall  FEL  extraction  (see  below).  This  can  be  overcome  by  altering  the 
resonance  condition.  Equation  {3.11},  along  the  length  of  the  undulator.  While  it  is 
possible  to  do  this  by  changing  the  undulator  period  \ ,  it  is  more  common  to  alter  K. 


This  typically  involves  changing  the  gap  between  undulator  poles  in  order  to  change  the 
magnetic  field  B,  so  the  process  is  usually  referred  to  as  “tapering”  the  undulator.  The 
effect  is  to  add  an  additional  phase  acceleration  term  to  the  pendulum  equation,  giving 


=  (?  -h  |a|  cos(i^  -I-  (p),  where 


S  =  -47rN- 


l  +  K^ 


AK 

K 


{3.28} 


The  quantity  S  is  referred  to  as  the  taper  rate[4].  For  taper  rates  that  are  too  large 
(<J>|a|),  the  electrons  are  accelerated  away  from  resonance  before  they  are  able  to 

bunch,  limiting  energy  transfer.  Rates  that  are  too  low  ( <?  <  4.^^ )  will  provide  only  a 

weak  effect,  and  do  not  enhance  the  FEL’s  natural  extraction.  Since  oscillators  typically 
have  undulators  too  short  to  face  this  problem,  tapering  is  much  more  commonly  applied 
to  FEL  amplifiers  [5]. 


D.  EXTRACTION 

Up  to  this  point,  references  to  “efficiency”  have  applied  to  the  overall  system.  To 
examine  specifically  the  efficiency  of  the  energy  exchange  between  the  electron  beam 
and  optical  field,  one  may  define  the  “extraction”  as  the  relative  amount  of  energy 
transferred  from  the  electron  beam  to  the  optical  field  on  a  single  pass,  that  is. 
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{3.29} 


E,  E,  y  ' 


where  is  the  energy  in  the  optical  field,  E^  is  that  in  the  electron  beam,  and 


Afi'  =  Aymc^  is  the  change  in  the  electron  beam  energy  over  the  length  of  the  undulator. 
Since  the  theoretical  PEL  gain  curve  has  a  width  of  A  v  =  2;r ,  it  follows  from  Equation 
{3.13}  that 


(A7)_(Av)^^  1 
y  AttN  2N 


{3.30} 


Equation  {3.30},  however,  applies  only  to  untapered  undulators. 
changes  the  expression  to  [5] 

1^1 

S;rN  ■ 


Introducing  taper 

{3.31} 


E.  BEAM  QUALITY 

The  derivations  of  the  FEE  pendulum  Equation,  {3.14},  and  the  FEE  wave 
Equation,  {3.27},  assumed  that  the  optical  fields  were  plane  waves  of  infinite  transverse 
extent.  This  assumption  remains  valid  mathematically  so  long  as  the  electron  beam 
remains  close  enough  to  the  center  of  the  optical  field  and  diffraction  is  small.  This 
places  a  condition  both  on  the  beam’s  radius  and  on  the  angular  spread  of  the  electrons. 
Neither  well-collimated  beams  of  large  radius  nor  highly  divergent  beams  of  small  initial 
radius  will  satisfy  the  assumptions  of  Parts  A  and  B.  One  appropriate  measure  of  beam 
quality,  then,  is  the  transverse  emittance,  defined  as 

£  =  {3.32} 

where  is  the  root-mean-square  (rms)  radius  of  the  electron  beam  and  is  its  rms 
angular  spread,  as  illustrated  in  Figure  4. 
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a)  Well-Collimated  Beam  b)  Highly  Divergent  Beam 

n 


Figure  4.  Beam  Quality  Dimensions 

An  alternate  measure,  which  has  the  advantage  of  being  conserved  as  the  beam  is 
accelerated,  is  the  normalized  transverse  emittance,  defined  as 

Sn=Yh0,.  {3.33} 

Equation  {3.13}  shows  that  changes  in  the  phase  velocity  are  proportional  to  the 
relative  change  in  Lorentz  factor,  dy !  y .  This  is  nonzero  for  a  relativistic  beam  from  a 
particle  accelerator  because  the  electrons  do  not  all  have  the  same  kinetic  energy,  and  is 
referred  to  as  the  “energy  spread.”  It  is  another  important  measure  of  the  beam’s  quality, 
and  it  influences  the  FEE  interaction,  as  well  as  the  beam  transport  devices  (e.g.,  focusing 
and  defocusing  elements). 
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IV.  SIMULATION  MODELS 


A.  FEL3D/FEL4D 

Faculty  and  students  at  the  Naval  Postgraduate  School  have  been  engaged  in 
theoretical  research  and  numerical  modeling  of  FELs  for  over  twenty  years.  Their  efforts 
have  supported  the  FEE  research  community’s  full  spectrum  of  activity — from  design 
support  for  new  experiments  to  refining  the  understanding  of  experimental  results.  Over 
time,  the  Naval  Postgraduate  School  (NPS)  has  developed  a  suite  of  simulations  for  both 
amplifiers  and  oscillators  that  have  been  benchmarked  against  both  theoretical  formulas 
and  experimental  results.  The  suite  includes  four  types  of  programs:  single-mode, 
longitudinal  multimode,  transverse  multimode  and  four-dimensional  [5]. 

The  single-mode  simulations  follow  the  phase-space  evolution  of  many  electrons 
and  one  optical  mode  (an  infinite  plane  wave)  along  the  undulator,  using  the  pendulum 
and  wave  Equations  ({3.14}  and  {3.27})  to  self-consistently  determine  the  behavior  of 
both  the  electrons  and  the  light.  They  can  be  used  to  independently  model  certain 
physical  effects,  allowing  for  general  understanding  of  the  origins  and  consequences  of 
those  effects.  They  are  small  and  computationally  inexpensive,  easy  to  understand,  and 
provide  a  solid  foundation  from  which  to  build  more  comprehensive  simulations. 

The  longitudinal  multimode  simulations  are  used  to  examine  the  evolution  of 
longitudinal  modes  in  the  optical  field.  The  optical  amplitude  and  phase  are  allowed  to 
vary  as  a  function  of  z  in  each  longitudinal  “slice”  of  the  optical  field.  Longitudinal 
modes  play  an  important  role  in  some  FELs  due  to  the  phenomenon  of  “pulse  slippage.” 
When  the  pulses  are  short,  the  difference  in  velocities  between  electron  and  optical  pulses 
can  significantly  affect  the  manner  in  which  the  light  is  amplified. 

The  transverse  multimode  simulations,  collectively  referred  to  as  FEL3D,  break 
from  the  assumption  that  the  optical  fields  are  plane  waves.  They  use  the  pendulum  and 
wave  equations  to  self-consistently  follow  the  evolution  of  the  electron  beam  and  the 
optical  field  in  x,  y,  and  r.  They  also  allow  for  axial  asymmetries,  such  as  separate 
emittances  in  the  two  transverse  planes,  astigmatism  in  oscillator  mirrors,  and  shifts  or 
tilts  of  either  the  electron  beam  or  the  oscillator  mirrors.  These  effects  can  also  have  a 
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significant  impact  on  PEL  gain  and  extraction.  Both  the  transverse  and  longitudinal 
multimode  programs  require  only  moderate  computing  resources,  and  run  on  personal 
computers,  though  their  efficiency  is  increased  when  multiple  processors  are  available. 

To  provide  a  more  realistic  simulation  by  following  both  longitudinal  and 
transverse  effects,  NPS  has  recently  developed  the  FEL4D  programs.  These  model  the 
electrons  and  the  light  in  x,  y,  and  z  while  allowing  them  to  evolve  in  time.  These 
simulations  require  significant  processing  time  and  large  amounts  of  memory,  and  are 
used  with  the  dedicated  resources  of  a  cluster  computer. 

B.  GINGER 

The  history  of  the  code  Ginger  began  with  the  code  Fred,  a  two-dimensional  PEL 
amplifier  simulation  code  developed  at  Lawrence  Livermore  National  Laboratory 
(LLNL)  some  three  decades  ago.  The  early  spinoffs  from  Fred  included  a  harmonic 
simulation  called  Nutmeg  and  a  three-dimensional  (i.e.,  multimode)  program  dubbed 
FredSD. 

Researchers  at  LLNL  and  Lawrence  Berkeley  National  Laboratory  (LBNL) 
created  Ginger  in  the  mid-1980s  primarily  to  examine  the  physics  and  consequences  of 
sideband  growth  in  single-pass  PEL  amplifiers  [6].  Ginger  is  polychromatic,  time- 
dependent  and  multidimensional.  It  has  grown  over  the  years  to  include  many  optional 
features,  such  as  a  shot-noise  startup  model,  configurations  for  complicated  PEL  types 
(e.g.,  harmonic  cascade  PELs),  and  higher-harmonic  emission  models.  Historically,  it 
has  been  limited  to  axisymmetric  (r,  z,  f)  geometry  for  fields,  but  recent  modifications 
include  optional  “slab”  geometry  (y,  z,  t,  infinite  x).  The  simulation  can  also  still  be  run 
in  a  non-time-dependent  (3D)  configuration,  referred  to  as  “FreJ  mode.” 

Ginger  and  its  postprocessor  are  both  flexible  and  powerful.  Both  are  written  in 
Fortran90  and  originally  designed  for  use  with  Unix,  but  have  been  successfully  ported 
to  various  platforms.  The  simulations  explored  for  this  thesis  operate  efficiently  on  a 
desktop  computer;  larger  and  more  complicated  problems  can  take  advantage  of  the 
multiprocessor  version  of  the  program  for  use  on  cluster  computers. 
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C.  FEL3D  CASE  STUDY:  BEAM  QUALITY  IN  A  GENERIC  EEL 
OSCILLATOR 

Table  1  represents  a  parameter  set  for  a  generie  PEL  oscillator.  These  values  do 
not  correspond  to  any  given  machine,  but  are  generally  illustrative  of  realistic  values  and, 


thus,  useful  as  a  case  with  which  to  explore  a  simulation  code’s  operation. 


Electron  Beam  Parameters 

Energy: 

Eb=\2Q  MeV 

Radius: 

=  0.1  mm 

Bunch  charge: 

^  =  1  nC 

Bunch  duration: 

tb  =2  ps 

Trans,  emittance: 

s„=  10  mm-mrad 

Long,  emittance: 

Si  =  200  keV-ps 

Peak  current: 

Ipeak=500A 

Energy  Spread: 

Ay/y  =0.2% 

Undulator  Parameters 

Period: 

To  =  3  cm 

Length: 

L  =  0.75  m 

Gap 

g=lcm 

Number  of  periods: 

N=25 

Undulator  parameter: 

^rms=1.78 

Resonator  Parameters 

Wavelength: 

A=  1.12  pm 

Length: 

S=20m 

Rayleigh  length: 

Zo  =  7.5  cm 

Quality  factor: 

Qn  =  4 

Table  1.  Generic  PEL  Oscillator  Parameters 


Most  of  these  quantities  are  self-explanatory  or  have  been  discussed  before. 
Longitudinal  emittance  is  an  electron  beam  quality  measure  that  is  proportional  to  energy 
spread.  The  undulator  gap  g  is  the  spacing  between  magnet  poles,  determining  the 
magnetic  field  strength  and,  thus,  the  undulator  parameter  K.  The  Rayleigh  length  Zq  is 
the  distance  from  the  optical  beam  waist  over  which  the  light  will  propagate  before  its 
cross-sectional  area  doubles.  The  optical  wavelength  value  shown  is  the  resonant 
wavelength  of  the  PEL,  determined  from  beam  energy  and  undulator  parameters  by 
Equation  {3.11}  (recall  that  =  ymc^).  The  quality  factor  Q„  represents  the  losses  per 

pass  in  the  optical  resonator,  with  low  Q„  values  corresponding  to  high  losses.  In  a 
conventional  laser,  this  includes  in-cavity  losses  such  as  heating  in  the  lasing  medium. 
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For  an  FEL  oscillator,  it  is  used  to  indicate  the  amount  of  power  outcoupled  through  the 
resonator’s  partially  transmissive  mirror.  Ideally,  it  is  related  to  the  outcoupling  mirror’s 
refleetivity  R  by 

R  =  l-—  {4.1} 

a 

The  output  of  FEL3D  is  mostly  graphieal,  though  speeifie  numeric  quantities  ean 
be  extracted  with  simple  Unix  commands.  Figure  5  is  part  of  the  output  for  the  ease 
outlined  in  Table  1,  depleting  the  x  =  0sliee  of  the  optical  field  after  300  passes.  The 
optical  field  amplitude  a(0,y,r)  is  shown,  with  lighter  blue  meaning  more  intensity. 

The  thin  cyan  lines  correspond  to  the  e"'  points  relative  to  the  peak  amplitude.  Spaeed 
along  the  x  axis  are  sample  eleetrons  (red  dots) — the  faet  that  they  all  fall  within  the  inner 
portion  of  the  optical  field  indieates  that  the  beam  and  field  were  well-matehed.  Also 
depicted  at  r  =  0  and  r  =  1  are  line-plots  of  the  optical  mode  shape,  with  magenta 
corresponding  to  the  eavity  fundamental  mode  for  referenee  and  yellow  to  the  mode 
observed  in  the  simulation.  Some  deviation  between  the  two  is  typical,  especially  at 
r  =  1 ,  and  the  “hollow”  (depressed-peak)  mode  shown  here  is  not  uncommon. 


Figure  5.  Generic  Oscillator:  Centerline  Optieal  Field  and  Electron  Beam  (FEL3D) 


Figure  6.  Generie  Ocillator:  Power  and  Gain  Evolution  (EEL 3D) 
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Another  portion  of  the  output  is  displayed  in  Figure  6.  This  depicts  the  evolution 
of  outcoupled  optical  power  and  optical  gain  over  300  passes,  in  normalized  units.  Both 
quantities  exhibit  transient  fluctuations  for  the  first  100  or  so  passes,  then  reach  steady- 
state. 

Several  other  figures  are  produced  for  a  given  case  in  FEL3D  and  can  be 
reviewed  either  to  check  simulation  viability  or  to  look  for  particular  behaviors.  Often, 
however,  the  goal  of  a  specific  design  (i.e.,  parameter  set)  is  a  particular  quantity,  usually 
the  output  power  or  extraction.  To  demonstrate  the  FEL3D  code  and  gain  a  general 
insight  into  the  behavior  of  FEL  oscillators,  we  conducted  a  parameter  study  to  determine 
the  extraction  against  simultaneous  variation  of  transverse  emittance  and  energy  spread. 
All  other  parameters  used  the  values  in  Table  1. 


Figure  7.  Generic  Oscillator:  Beam  Quality  Parameter  Study  (FEL3D) 

The  results  of  this  study  are  shown  in  Figure  7.  The  results  point  to  the  general 
conclusion  that  increasing  either  emittance  or  energy  spread  reduces  extraction.  This  is 
somewhat  intuitive,  but  studies  such  as  these  can  quantify  the  dependence.  For  these 
parameters,  the  peak  extraction  is  approximately  3%,  a  strong  but  plausible  value.  The 
two  red  contours  on  this  graph  correspond  to  2%  and  1%  extraction.  Staying  above  half 
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the  peak  extraction  requires  simultaneously  keeping  the  emittance  under  30mm-mrad  and 
the  energy  spread  under  0.4%.  Both  requirements  are  reasonable  values  by  current 
standards. 

D.  GINGER  CASE  STUDY:  JLAB’S  14KW  EXPERIMENT 

In  October  2006,  Thomas  Jefferson  National  Laboratory  (JLab)  broke  their  record 
for  PEL  average  output  power  by  achieving  14.2kW  in  continuous-wave  oscillator 


operation.  The  parameters  for  that  experiment  are  summarized  in  Table  2. 


1  Electron  Beam  Parameters  | 

Energy: 

Eb=  150  MeV 

Radius: 

rb  =  0.17  mm 

Bunch  charge: 

114pC 

Bunch  duration: 

tb  =0.4  ps 

Peak  current: 

Ipeak  =  285  A 

Avg  current: 

^ave  8.5inA 

Trans,  emittance: 

£•„  =  8  mm-mrad 

Long,  emittance: 

Si  =  80  keV-ps 

Energy  Spread: 

Ay/y  =0.2% 

Undulator  Parameters 

Period: 

Ao  =  5.5  cm 

Length: 

Z=  L65m 

Number  of  periods: 

N=30 

Undulator  parameter: 

Krms=l.36 

Resonator  Parameters 

Length: 

S=32m 

Wavelength: 

A=  1.6  pm 

Rayleigh  length: 

Zo  =  75  cm 

Waist  radius: 

Wo  =  0.6  mm 

Mode  radius  at  mirrors: 

w  =  1.3  cm 

Quality  factor: 

Qn  =  5 

Table  2.  JLab  14kW  Parameters 


These  parameters  were  used  to  build  familiarity  with  Ginger  just  as  the  generic 
PEL  oscillator  parameters  were  used  to  build  familiarity  with  FEL3D.  A  brief 
comparison  shows  that  the  numbers  in  Tables  1  and  2  are  of  similar  magnitude.  More 
importantly,  Table  2  represents  a  known  experiment  with  known  results,  which 
furthermore  has  been  previously  simulated  by  NPS. 
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The  output  of  Ginger's  dedicated  postprocessor  is  a  series  of  graphs.  Most  of 
them  track  different  quantities  than  the  output  of  FEL3D,  reflecting  the  difference  in  the 
programs’  origins.  Gain  and  power  evolution,  however,  are  important  to  any  PEL  and 
appear  in  both  codes’  output.  For  the  case  outlined  in  Table  2,  the  gain  evolution  is 
shown  in  Figure  8. 


Total  Gain  vs.  PASS  # 

1.0 

0.6 
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1  41  81  120  160  200 

PASS# 

JLab  14kW:  Gain  Evolution  (Ginger) 

As  with  the  previous  example,  the  simulation  shows  significant  transients  over  the 
first  one-third  of  the  evolution,  then  settles  into  steady-state.  The  steady-state  gain  value 
appears  to  be  approximately  25%.  The  gain  per  pass  G  required  for  steady-state 
operation  is  given  by 

G  =  ^— ,  {4.2} 

a-1 

or  25%  for  Q„=5 .  This  is  a  simple  check  that  the  simulation  is  behaving  as  expected 
and  increases  confidence  in  its  results. 

The  power  evolution  is  depicted  in  Figure  9.  The  transients  are  minimal — a 
simple  ramp-up — and  the  FEE  reaches  steady-state  around  pass  60,  just  as  the  gain  did. 
This  is  another  check  on  behavior  meeting  expectations.  The  values  depicted  on  this 
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Figure  8. 
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graph  are  of  peak  power  at  the  exit  of  the  undulator  To  eonvert  to  average 

outeoupled  power  one  must  know  the  pulse  duration  tb,  the  pulse  repetition 

frequeney  Q  and  quality  faetor  Q„.  The  formula  relating  them  is 


P  =P 

out  und 


Q, 


{4.3} 


V  J 


For  this  experiment,  the  quality  faetor  and  pulse  duration  are  given  in  Table  2,  while  the 
PRF  is  approximately  75MHz.  Using  1.7GW  for  P^^^  yields  P^^^  « 10.2kW  . 


Output  Radiation  Power  vs.  PASS  # 

2.0E+09 

1.6E+09 

^1.2E+09 

ts 

0) 

^8.0E+08 
O 
Q. 

4.0E+08 

O.OE+00 

1  41  81  120  160  200 

PASS# 

Figure  9.  JLab  14kW:  Power  Evolution  {Ginger) 

While  this  result  is  approximately  27%  below  the  known  experimental  results,  it 
is  of  the  eorreet  order  of  magnitude.  Furthermore,  previously  eondueted  NFS  3D 
simulations  predieted  nearly  20kW,  a  40%  overestimation.  The  eonelusion  drawn  at  the 
time  was  that  pulse  slippage  effeets  were  present  in  the  experiment,  whieh  would  reduee 
the  predieted  output  power. 
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E.  COMPARATIVE  CASE  STUDY:  BEAM  QUALITY  IN  A  GENERIC  EEL 

AMPLIFIER 

To  compare  the  results  of  FEL3D  and  Ginger,  we  took  the  ease  of  a  generic  FEL 
amplifier,  using  the  parameters  in  Table  3.  The  values  ehosen  mateh  closely  with  the 
generic  FEL  oscillator  parameters  from  Part  C  by  design.  The  most  notable  difference  is 
the  significant  increase  in  undulator  length,  typieal  for  amplifier  designs.  Additionally, 
the  optieal  eavity  parameters  have  been  replaced  by  seed  laser  speeifieations.  The  seed 
wavelength,  whieh  will  also  be  the  output  wavelength,  is  determined  entirely  by  the  type 
of  seed  laser  used  and  therefore  typically  driven  by  the  needs  of  the  application.  In  this 


case,  the  value  was  ehosen  to  match  the  output  wavelength  of  the  generie  oseillator. 


1  Electron  Beam  Parameters  | 

Energy: 

Eb  =120  MeV 

Radius: 

Tb  =0.3  mm 

Buneh  charge: 

q=l  nC 

Buneh  Duration: 

tb  =2  ps 

Trans.  Emit: 

8n  =  10  mm-mrad 

Long,  emit: 

8i  =  200  keV-ps 

Peak  Current: 

Ipeak  =  500  A 

Energy  Spread: 

Ay/y  =0.2% 

Undulator  Parameters 

Period: 

A-o  =  3  em 

Length: 

L  =  3.6  m 

Gap 

1  cm 

Number  of  periods: 

N=  120 

Undulator  parameter: 

Krms=  1.78 

Seed  Laser  Parameters 

Wavelength: 

A  =  1.12  pm 

Pulse  Duration: 

topt  =  4  ps 

Average  power: 

Popt=  1  W 

Peak  power: 

Ppeak=330kW 

Table  3.  Generic  FEL  Amplifier  Parameters 


Using  FEL3D,  we  eondueted  a  parameter  study  on  beam  quality.  As  before,  we 
varied  transverse  emittanee  and  energy  spread  simultaneously  and  noted  output 
extraction.  For  this  version  of  the  study,  we  included  a  taper  rate  of  AB  /  B  =  5%  over 
the  last  meter  of  the  undulator.  The  results  are  shown  in  Figure  10. 
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Figure  10.  Generic  Amplifier:  Beam  Quality  Parameter  Study  (FEXiZ)) 

Peak  extraction  for  this  case  was  somewhat  lower  than  for  the  oscillator,  near  two 
percent.  This  is  a  normal  and  desirable  value  by  current  design  standards.  To  stay  above 
half  this  peak  value,  it  was  observed  that  the  design  must  maintain  transverse  emittance 
below  lOmm-mrad  and  energy  spread  below  0.25%.  Recall  that  the  corresponding  limits 
for  the  oscillator  were  30mm-mrad  and  0.4%.  This  amplifier,  then,  is  shown  to  be  more 
sensitive  to  beam  quality  than  a  similar  oscillator.  This  result  is  consistent  with  known 
experiments  and  the  underlying  physics.  Put  simply,  amplifiers  have  significantly  longer 
undulators;  therefore,  any  degradation  the  beam  quality  introduces  to  the  interaction  is 
compounded  over  a  longer  time  and  has  a  greater  effect.  In  an  oscillator,  each  pass  an 
optical  pulse  makes  through  the  resonator  is  a  new  interaction  with  new  electrons,  so 
beam  quality  effects  are  essentially  “reset”  each  time. 

For  the  code  comparison,  we  examined  the  variation  of  emittance  and  energy 
spread  separately.  Figure  1 1  shows  the  results  of  varying  emittance  while  holding  energy 
spread  constant  at  0.2%. 
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Extraction  vs.  Emittance 


8^  (mm-mrad) 


Figure  11.  Generic  Amplifier:  Extraction  vs.  Emittance  (Comparative) 

The  agreement  between  the  codes  is  close.  The  only  notable  discrepancy  occurs  for  the 
minimum  emittance  value.  Neither  simulation  included  taper  in  this  case;  the  inset  figure 
displays  the  power  evolution  from  FEL3D,  which  shows  clearly  that  the  interaction 
reached  saturation  before  the  end  of  the  undulator. 

Figure  12  shows  the  results  of  varying  energy  spread  while  holding  emittance  at 
lOmm-mrad.  Again,  the  codes  agree  very  well  until  low  values  of  energy  spread.  The 
inset  indicates  the  same  reasons  as  above  for  this  discrepancy. 


Extraction  vs.  Energy  Spread 


Ay/y  (%) 


Figure  12.  Generic  Amplifier:  Extraction  vs.  Energy  Spread  (Comparative) 
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F.  CONCLUSIONS 

The  two  simulation  codes  considered — FEL3D  and  Ginger  in  its  3D  mode — 
appear  consistent.  Some  variation  in  the  detailed  results  should  be  expected  due  to 
internal  differences  in  the  simulations.  The  strongest  disparity  seen  in  our  examination 
came  from  the  Jefferson  Lab  14kW  results,  which  Ginger  underestimated  and  FEL3D 
overestimated.  Both  results,  however,  were  of  the  correct  order,  and  both  discrepancies 
are  likely  due  to  the  lack  of  pulse-slippage  effects  in  the  models.  In  direct  comparison, 
the  two  codes  agreed  very  well  on  beam-quality  effects  in  a  generic  FEL  amplifier.  The 
discrepancies  seen  here  show  that  the  NFS  models  reached  saturation  early  and  therefore 
exhibited  a  lower  extraction  than  they  would  have  with  a  taper. 
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V.  ELECTROSTATIC  CATHODE  TEST  CELL  DESIGN 


As  discussed  in  Chapter  II,  design  of  an  injeetor  is  important  to  overall  design  of  a 
free  eleetron  laser  (PEL).  Characterization  of  the  cathode  used  is  important  to  the  design 
of  any  injector.  Normally,  since  injectors  are  generally  not  built  to  support  more  than  one 
cathode  type,  this  is  performed  on  a  “one-off’  basis  for  the  eathodes  used  by  a  given 
PEL.  Additionally,  replaeing  a  eathode  is  often  an  involved  task,  undertaken  only  when 
the  one  in  use  has  been  spent  or  damaged. 

A  tool  the  PEL  community  finds  useful  is  the  eathode  test  eell:  a  deviee  that  is 
not  a  full-fledged  injeetor  but  can  be  used  to  test  and  eharacterize  different  cathodes 
under  quasi-operational  eonditions.  Due  to  the  signifieant  differenees  between  the 
physical  nature  of  and  measurements  required  for  thermionie,  photoeleetric,  and  field 
emission  eathodes,  the  test  eell  design  presented  in  this  thesis  foeuses  only  on  one  type, 
namely  field  emission  (PE).  Existing  PE  test  cells  inelude  the  one  used  to  explore 
diamond  field-emitter  array  (DPEA)  eathodes  at  Vanderbilt  University  (figure  13). 
Typically,  test  eells  aehieve  high  fields  with  low  voltages  using  small  anode-cathode 
spaeing  (“gap”).  At  higher  voltages,  the  small  gap  and  gridded  anode  ean  result  in 
subjeeting  the  device  and  the  eathode  to  significant  damage  from  vaeuum  arcing  and  ion 
blowback. 


Collector: 


figure  13.  Vanderbilt’s  DFEA  Test  Stand  [From  7] 
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The  general  goal  for  the  design  presented  in  this  chapter  is  for  a  higher-voltage, 
larger-gap,  non-gridded  test  cell.  Building  such  a  device  will  allow  NFS  researchers  to 
explore  the  properties  of  FE  cathodes  at  higher  voltages  and  gain  general  experience  with 
FE  cathodes  for  possible  use  in  their  future  FEE  experimental  facility.  The  design  work 
discussed  here  presents  preliminary  considerations — building  a  foundation  upon  which 
future  research  will  build,  until  a  final  design  is  selected  and  the  system  procured. 

A.  FIELD  EMISSION  THEORY 

Cold  (i.e.,  not  thermionic)  field  electron  emission  can  occur  at  the  surface  of  a 
bulk  material  when  strong  static  external  electric  fields  are  applied.  It  arises  from  the 
phenomenon  of  quantum  tunneling  and  its  explanation  in  the  1920s  was  one  of  the  early 
triumphs  in  the  then-emerging  field  of  quantum  mechanics.  (It  is  also  important  due  to 
its  many  practical  applications,  included  devices  such  as  the  test  cell  under  discussion.) 
The  current  density  J  emitted  in  the  presence  of  applied  field  E  is  given  by  the  Fowler- 
Nordheim  equation,  a  simplified  form  of  which  is  [8] 

J  =  AEV'''^E  {5.1} 

The  constants  A  (in  units  of  A/V^)  and  B  (in  units  of  E  /  m )  are  related  to  the  material’s 
work  function  (p  and  Fermi  energy  by 

^  _  6.2-1 0~^ 

(p  +  Sj  ]l  ^  {5.2} 

5  =  (6.8-10^)^3^ 

Several  other,  more  complicated,  forms  of  the  Fowler-Nordheim  equation  exist,  since 
several  methods  have  been  used  to  derive  the  fundamental  form  of  the  coefficients  A  and 

B.  Equation  {5.1}  is  presented  in  this  form  to  illustrate  the  most  salient  points  of  the 
formula:  inverse  exponential  dependence  on  E  and  variation  with  material  properties. 
The  fundamental  principle  is  that  a  larger  field  will  yield  more  current. 

When  the  surface  of  the  bulk  material  is  not  perfectly  planar,  any  protrusions  will 
cause  local  field  enhancements  (field  lines  will  converge  at  their  tips),  effectively 
multiplying  the  applied  field  by  a  geometry-dependent  factor  at  those  points.  This 
significantly  lowers  the  overall  field  E  required  to  reach  a  given  value  of  current  density 
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J.  The  field  E  as  used  in  Equation  {5.1}  remains  the  pertinent  value  at  the  tips.  This 
occurs  naturally  due  to  the  microprotrusions  (i.e.,  surface  roughness)  of  bulk  materials, 
but  it  can  also  be  deliberately  exploited  by  building  in  macroscopic  protrusions. 

Diamond  field  emitter  arrays  (see  Figure  14)  take  advantage  of  this  technique.  A 
thin  layer  of  diamond  is  grown  on  a  wafer  of  a  suitable  substrate  such  that  the  diamond 
forms  into  an  array  of  pyramids.  When  placed  in  an  external  field,  the  pyramids  serve  as 
field-enhancement  points  and  can  emit  electrons,  which  will  be  accelerated  by  the  applied 
field  and  form  beamlets.  Arrays  can  be  fairly  dense  (in  terms  of  tips  per  unit  area)  to 
provide  high  current  density,  but  this  will  cause  space-charge  forces  to  affect  the 
beamlets’  propagation. 


Figure  14.  Example  of  a  DFEA  [After  9] 

B.  TEST  CELL  DESIGN  REQUIREMENTS 

A  cathode  test  cell  must  apply  a  given  field  to  the  cathode  being  tested  and  allows 
characterization  or  monitoring  of  the  cathode’s  performance.  Proper  characterization  of 
an  FE  cathode  involves  measuring  the  current  emitted  as  a  function  of  field  and  time,  and 
monitoring  the  beam  shape,  which  will  have  an  effect  on  beam  transport  and  emittance. 
For  DFEA  cathodes,  this  is  sometimes  analyzed  on  a  per-tip  basis;  for  sparse  enough 
arrays,  the  beam  from  each  separate  tip  can  be  monitored.  To  support  testing  of  multiple 
cathodes,  the  cathode  holder  should  be  easily  removable.  To  facilitate  varying  the  field 
strength  E,  either  the  power  supply  voltage  V  must  be  variable,  or  the  cathode  assembly 
should  be  movable,  allowing  for  adjustment  of  the  anode-cathode  gap  d,  or  both.  For  the 
ideal  case — where  the  anode  and  cathode  are  infinite  parallel  plates — the  three  quantities 
are  related  by 
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E  =  Vld.  {5.3} 

Operation  of  a  field-emission  (FE)  cathode  requires  a  vacuum  chamber — any  gas 
present  would  be  ionized  by  the  electron  beam  and  likely  damage  the  cathode.  Figure 
15  displays  the  insulating  ceramic  upon  which  this  test  cell  design  was  based.  Its  metal 
flanges  are  separated  by  a  ceramic  break,  allowing  the  components  at  either  end  to  be 
electrically  isolated  from  one  another.  The  flanges  are  of  a  standard  type  (“conflat”), 
compatible  with  many  common  devices  and  fittings.  Custom  fittings  can  easily  be 
ordered  once  the  required  configuration  of  measurement  devices,  power  ports,  vacuum 
system  ports,  etc.,  are  fully  determined. 


Figure  15.  Ceramic  Insulator  with  Conflat  Flanges 

The  electron  beam’s  size  can  be  analyzed  by  placing  a  phosphor  screen  in  the 
beam’s  path  and  observing  the  spot(s)  illuminated  upon  it.  The  phosphor  screen  can  also 
serve  as  a  charge  collector,  with  current  easily  measured  across  a  resistor  placed  between 
it  and  ground.  The  stability  of  the  current  measured  can  be  monitored  using  a  current 
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transformer  placed  around  any  isolated  portion  of  the  chamber  through  which  the  beam 
travels.  Using  the  type  of  cylinder  depicted  in  Figure  15,  one  can  conceive  of  an 
arrangement  such  as  shown  in  Figure  16. 


Here,  the  movable  cathode  assembly  is  fixed  to  an  endcap  that  is  bolted  to  the 
lower  flange,  allowing  for  easy  removal.  This  end  is  held  at  a  large  negative  voltage 
during  operation.  The  annular  anode  is  electrically  coupled  to  the  other  endcap,  which 
contains  the  phosphor  screen  (beam  target)  and  is  grounded  during  operation.  Figure  16 
provides  the  basis  for  the  preliminary  test  cell  design  presented  in  the  following  sections. 
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Figure  16. 


Electrostatic  Cathode  Test  Cell  Schematic 
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C.  DESIGN  TRADEOFFS  AND  PARAMETER  SELECTION 

Experimentally,  the  emission  predicted  by  Equation  {5.1}  is  generally  observable 
only  once  the  applied  field  E  has  been  raised  to  some  minimum  threshold  value, 
dependent  on  the  geometry  and  material.  For  the  purposes  of  this  thesis,  the  threshold 
field  is  considered  to  be  lOMV/m,  and  the  maximum  operating  field  is  20MV/m.  These 
numbers  are  roughly  characteristic  of  the  DFEA  cathodes  to  be  studied  at  NFS. 

In  addition  to  the  insulating  ceramic  pictured  above,  NFS  has  on  hand  a  fixed- 
voltage  lOOkV  power  supply  suitable  for  use  with  this  test  stand.  Field  strength 
adjustment  will  therefore  be  made  through  adjustment  of  the  anode-cathode  gap.  The 
pertinent  tradeoff  is  between  field  strength  and  gap  size:  operating  with  a  reduced  the 
gap  will  achieve  higher  field  but  risk  causing  vacuum  arcs,  potentially  damaging  the  test 
cell  and  cathode. 

This  tradeoff  also  affects  the  placement  of  the  phosphor  screen.  The  simplest 
arrangement  would  be  to  fix  the  screen  directly  to  the  anode,  but  this  creates  significant 
risk  of  ion  blowback  as  it  is  bombarded:  ions  excited  off  the  screen  would  be  positive 
and  therefore  be  accelerated  toward  the  cathode.  This  risk  can  be  reduced  by  placing  the 
screen  past  the  anode  in  an  approximately  field-free  region;  any  ions  liberated  from  the 
screen  will  remain  in  its  vicinity  due  to  the  lack  of  field  instead  of  accelerating  back  to 
damage  the  cathode.  This  arrangement,  however,  requires  an  aperture  in  the  anode, 
which  will  alter  the  field  distribution  and  potentially  defocus  the  beam.  Some  setups 
(including  that  in  Figure  13)  mitigate  this  by  spanning  the  aperture  with  a  grid  or 
“pepperpot,”  which  mostly  allows  the  electrons  to  pass,  but  keeps  the  field  roughly 
uniform.  This  increases  the  risk  of  damaging  the  cathode  and  anode,  however,  albeit  less 
so  than  with  the  direct- impingement  arrangement.  The  grid’s  focusing  effects  can  also 
lead  to  emittance  degradation.  For  this  design,  the  aperture  will  remain  ungridded; 
avoiding  arcing  and  ion  blowback  is  assumed  to  outweigh  the  effects  of  defocusing.  This 
tradeoff  remains  a  viable  field  of  study  for  future  design  work. 

Though  the  electron  beam  will  emerge  from  the  cathode  approximately  collimated 
and  perpendicular  to  the  cathode  surface,  space-charge  forces  will  tend  to  defocus  the 
beam  as  it  propagates.  Techniques  exist  to  counteract  this,  including  the  famous  “Fierce 
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geometry”  [10].  In  a  Pierce  configuration,  the  metal  surface  on  which  the  cathode  rests  is 
angled  in  the  direction  of  beam  travel,  forcing  the  local  field  lines  to  point  inward  and 
help  focus  the  beam.  This  is  shown  in  Figure  17,  where  purple  lines  represent 
equipotentials  and  red  arrows  are  local  field  vectors.  The  Pierce  effect  can  help  mitigate 
space-charge  effects,  but  since  it  inherently  enhances  the  field  in  certain  regions,  it 
increases  the  risk  of  breakdown.  Pierce  geometries  also  reduce  the  field  at  the  cathode 
surface  however  (as  is  apparent  in  Figure  17),  requiring  a  smaller  gap  for  the  same 
voltage.  In  this  design,  Pierce  geometry  is  included  at  the  cathode  support  but  not 
analyzed  in  detail.  The  topic  remains  open  for  optimization  in  future  studies. 


Figure  17.  An  Example  of  Pierce  Geometry  [After  1 1] 


The  design’s  major  parameters  are  summarized  in  Table  4.  The  voltage  value  is 
fixed  by  available  equipment,  and  the  field  strengths  represent  design  goals;  the  range  of 
gap  values  is  the  subject  of  study  below,  with  initial  estimates  based  on  Equation  {5.3}. 


Parameter 

Value 

Operating  Voltage 

lOOkV 

Anode-Cathode  Gap 

5-2. 5mm 

Operating  Field 

10-20MV/m 

Table  4. 

Guiding  Design  Parameters 
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D. 


THE  DESIGN  TOOL:  POISSON  SUPERFISH 


First  developed  at  Los  Alamos  National  Laboratory  in  the  late  1970s,  Poisson 
Superfish  is  a  suite  of  programs  used  to  calculate  fields  (electrostatic,  magnetostatic  or 
RF  electromagnetic)  in  two-dimensional  geometries.  The  component  programs  used  for 
this  design  work  include  Automesh,  which  generates  a  triangular  mesh  (of  variable 
spacing)  fitted  to  the  boundaries  of  the  different  regions  of  the  geometry;  Poisson,  which 
calculates  static  fields  by  successive  over-relaxation;  the  WSFPlot  display  tool,  used  to 
generate  the  diagrams  shown  below;  and  SF7,  a  postprocessor  used  to  extract  specific 
information  from  the  field  map  (e.g.,  field  strength  along  a  specified  line). 

Figure  18  displays  an  example  of  WSFPlot' s  rendering  of  the  region  of  the 
schematic  indicated.  Input  geometries  are  specified  in  Cartesian  coordinates,  mapping 
x^r  and  y^z  for  cylindrically  symmetric  problems  such  as  this.  This  is  a  common 
approach,  and  a  flag  within  the  input  file  will  cause  output  tables  to  be  presented  as  (r,z) 
instead  of  (x,y).  Output  files  consist  of  a  table  representing  the  field  map.  WSFplot  is 
used  to  overlay  this  on  the  problem  geometry,  and  includes  several  useful  display 
options,  such  as  field  arrows,  equipotential  contours  and  the  mesh  used. 


Figure  18.  Superfish  Rendering  of  Test  Cell  Interior 
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The  Superfish  suite  is  a  powerful  tool  for  a  wide  variety  of  applications,  and 
includes  several  additional  programs  not  applicable  to  this  design  effort.  There  are  also 
multiple  postprocessors,  including  tabular  display  tools.  For  this  thesis,  however,  all 
plots  were  generated  in  Matlab  using  data  from  SF7. 

E.  ANODE  SHAPE 

The  baseline  concept  for  an  anode  in  this  design  is  a  simple  ring.  A  torus  of  metal 
with  rectangular  cross-section  would  be  easily  manufactured  and  serve  the  essential 
purposes  of  an  anode.  Figure  19  shows  an  initial  design  concept  using  this  type  of  anode 
as  a  close-up  of  the  anode-cathode  region.  It  is  a  portion  of  a  cross-sectional  drawing 
such  as  displayed  in  Figure  18  above;  the  left  border  is  the  test  cell’s  centerline.  The  red 
arrows  depict  the  magnitude  and  direction  of  the  field  at  given  points,  while  the  pink  lines 
are  equipotential  contours. 


The  anode  pictured  is  rectangular  in  cross-section,  with  its  lower  comer  is  filleted 
to  avoid  creating  a  field  enhancement.  This  principle  is  also  applied  to  the  outer  edge  of 
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the  cathode  support  stalk.  As  the  field  map  shows,  the  field  is  elevated  in  that  region 
even  with  the  fillet;  peak  field  value  there  should  be  monitored  in  future  design  iterations. 

The  cathode  is  pictured  in  black,  resting  inside  a  divot  in  order  to  have  Pierce 
focusing.  The  divot’s  comers  are  not  filleted,  but  do  not  appear  to  greatly  perturb  the 
field  distribution.  As  expected,  the  field  begins  to  defocus  near  the  anode  aperture,  as 
evidenced  by  the  sloping  equipotentials. 

Figure  20  shows  the  second  general  design  configuration  considered,  featuring  an 
angled  anode.  Angling  the  lower  surface  should  help  alleviate  the  field  concentration 
near  the  comer  of  the  cathode  support  stalk.  Angling  the  top  surface  helps  mitigate 
damage  to  the  anode  by  progressively  increasing  the  distance  between  it  and  the  electron 
beam.  If  the  anode’s  angle  is  greater  than  that  of  the  diverging  electron  beam,  no 
electrons  will  strike  (and  damage)  the  anode.  Both  of  these  features  will  affect  field 
distribution;  it  is  pertinent  to  ask  whether  they  measurably  affect  field  strength  at  the 
cathode. 
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F. 


GAP  AND  APERTURE  SPACING 


Given  a  fixed  voltage  for  the  power  supply,  we  studied  two  ways  of  altering  the 
field  distribution  at  the  cathode:  anode-cathode  gap  and  anode  aperture.  The  cathode’s 
own  geometry  has  an  effect  as  well,  but  was  not  studied  as  a  part  of  this  work.  As  a 
baseline,  the  initial  configuration  considered  for  each  anode  shape  uses  a  1cm  gap,  as 
measured  from  the  top  of  the  cathode  support  stalk  to  the  lowest  point  on  the  anode,  and  a 
1cm  aperture.  With  a  potential  difference  of  lOOkV,  this  gap  value  corresponds  to  a 
lOMV/m  gradient  for  infinite  parallel  plates  (see  Equation  {5.3}).  For  the  actual 
geometry,  the  results  are  significantly  lower,  as  shown  in  Figure  21,  due  to  Pierce  effects 
and  the  anode  aperture.  Field  magnitude  values  are  shown  as  a  function  of  radius  out  to 
5mm,  the  radius  of  the  cathode  divot.  The  cathode  itself  will  have  a  radius  of 
approximately  4mm.  This  baseline  configuration  data  exhibits  not  only  a  failure  to  reach 
the  minimum  operating  field  but  also  a  strong  radial  dependence  in  field  strength.  The 
fields  are  reasonably  uniform  only  for  the  first  l-2mm,  corresponding  to  25%  of  cathode 
area.  At  a  radius  of  2mm,  the  axial  field  is  reduced  by  2.5%  from  its  centerline  value  for 
both  anode  shapes. 


Figure  21.  Baseline  Design 

The  next  case  examined  keeps  the  same  gap  value  but  narrows  the  aperture  by 
half  If  this  showed  improvement  in  field  values,  it  could  present  significant  advantages. 
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chiefly  by  flattening  out  the  field  in  the  aperture.  As  Figure  22  shows,  however,  the 
improvement  is  slight,  showing  a  peak  field  of  6.6MV/m  vs.  6.4MV/m,  and  the  strong 
radial  dependence  remains.  As  before,  the  simple  ring  anode  provides  somewhat  more 
favorable  fields,  but  in  this  case  the  disparity  is  a  bit  greater.  Neither  design  appears 
overly  promising,  and  the  narrowed  aperture  introduces  the  same  risks  as  a  pepperpot 
design.  A  narrow  aperture  can  even  block  a  portion  of  the  beam.  While  this  is 
sometimes  used  in  injectors  to  reduce  beam  emittance,  it  is  not  useful  for  measuring  the 
cathode’s  characteristics,  and  could  damage  the  anode. 


Figure  22.  Narrow-aperture  Design 

The  next  design  iteration  returns  to  the  1cm  aperture  and  halves  the  gap.  The 
results,  seen  in  Figure  23,  show  significant  improvement.  The  peak  field  in  this  case 
actually  exceeds  the  target  value  of  lOMV/m.  The  radial  dependence  of  the  field  is 
somewhat  lower  but  still  significant,  reaching  the  2.5%  reduction  level  around  2.5mm. 
The  disparity  between  the  two  anode  shapes  has  also  narrowed,  indicating  that  the  angled 
anode  may  be  a  viable  option  at  this  gap  spacing.  While  this  case  meets  the  minimum 
requirements,  it  is  worth  investigating  yet  smaller  gap  sizes. 

Decreasing  the  gap  to  2.5mm  produces  the  results  shown  in  Figure  24.  For  this 
case,  peak  fields  fall  well  within  the  operating  range,  and  the  angled  anode  performs 
almost  as  well  as  the  simple  ring.  That  result  is  favorable  since  such  a  narrow  gap  carries 
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Figure  23.  Intermediate-gap  Design 


increased  risk  of  arcing  and  the  angled  anode  may  be  vital  to  lowering  fields  outside  the 
cathode  region.  Even  more  important  a  result  is  the  field’s  uniformity.  The  field  strength 
only  begins  to  drop  off  after  about  3.5mm,  meaning  the  field  is  uniform  within  2.5%  over 
75%  of  the  cathode’s  area. 


Figure  24.  Narrow-gap  Design 
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G.  RECOMMENDATIONS 


The  conflat-cylinder  insulator  appears  entirely  feasible  for  testing  cathodes.  From 
a  field-strength  perspective,  a  5mm  gap  is  recommended  as  the  maximum  operating 
setting,  and  2.5mm  as  the  nominal  setting,  using  a  1cm  or  wider  aperture.  Future  design 
work  should  examine  the  possibility  of  breakdown  in  the  gap  and  whether  an  angled 
anode  -  perfectly  suitable  for  achieving  the  required  fields  at  the  cathode — would  provide 
mitigation. 
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